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Abstract 

The bottomonium states due to their varying binding energies dissolve at different temperatures 

■ and thus their nuclear modification factors and relative yields have potential to map the properties 

■ of Quark Gluon Plasma (QGP). We estimate the suppression of bottomonia states due to color 
r~| . screening in an expanding QGP of finite lifetime and size with the conditions relevant for PbPb 

collisions at LHC. The properties of T states and recent results on their dissociation temperatures 
have been used as ingredient in the study. The nuclear modification factors and the ratios of yields 
of T states are then obtained as a function of transverse momentum and centrality. We compare 
our theoretical calculations with the bottomonia yields measured with CMS in PbPb collisions at 
y^SNN = 2.76 TeV. The model calculations explain the data very well. 
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I. INTRODUCTION 



The heavy ion colhsions produce matter at extreme temperatures and densities where it is 
expected to be in the form of Quark Gluon Plasma (QGP), a phase in which the quarks and 
gluons can move far beyond the size of a nucleon making color degrees of freedom dominant 
in the medium. The experimental effort in to produce such matter started with low energy 
CERN accelerator SPS and evolved through voluminous results from heavy ion collision at 
Relativistic Heavy Ion Collider (RHIC) [1]. The recent results from Large Hadron Collider 
(LHC) experiments ^ are pointing towards formation of high temperature system in many 
ways similar to the matter produced at RHIC. One of the most important signal of QGP is 
the suppression of quarkonium states [3], both of the charmonium {J/ip, ip{2S), Xc, etc) and 
the bottomonium (T(IS') , T(2S'), Xbi etc) families. This is thought to be a direct effect of 
deconfinement, when the binding potential between the constituents of a quarkonium state, 
a heavy quark and its antiquark, is screened by the colour charges of the surrounding light 
quarks and gluons. The ATLAS and CMS experiments have carried out detailed quarkonia 
measurements in Pb+Pb collisions with the higher energy and luminosity available at the 
LHC. The ATLAS measurements |4] show suppression of inclusive J/ip with high transverse 
momenta pt in central PbPb collisions compared to peripheral collisions at a/s^vjv ~ 2. 76 
TeV. Similarly, CMS measured a steady and smooth decrease of suppression of prompt J/ip 
as a function of centrality with nuclear modification factor Raa remaining < 1 even in the 



peripheral bin js]. 

The melting temperature of the quarkonia states depends on their binding energy. The 
ground states, J/ip and T(IS') are expected to dissolve at significantly higher temperatures 
than the more loosely bound excited states. The T(2S') and T(3S') have smaller binding 
energies as compared to ground state T(IS') and hence are expected to dissolve at a lower 
temperature. With the 2011 Pb+Pb run the CMS published results on sequential suppres- 
sion of T(nS') states as a function of centrality {g] with enlarged statistics over their first 
measurement j?] where a suppression of the excited T states with respect to the ground state 
have been observed in PbPb collisions compared to pp collisions at v^tvat ~ 2.76 TeV. 

The quarkonia yields in heavy ion collisions are also modified due to non-QGP effects 



such as shadowing, an effect due to change of the parton distribution functions inside the 
nucleus, and dissociation due to nuclear or comover interaction [8]. Due to higher mass, the 
nuclear suppression is expected to be less for bottomonia over charmonia. If large number of 
heavy quarks are produced in initial heavy ion collisions at LHC energy this could even lead 
to enhancement of quarkonia via statistical recombination jo], [lo|. The effect of regeneration 
is expected to be less significant for bottmonia as compared to charmonia since bottom 
quarks are much smaller in number as compared to charm quarks. In addition, due to higher 
bottom mass the bound state properties obtained from potential models are more reliable. 
Thus recent years witness a shift in the interest to bottomonia. The ratios of the yields of 
excited states to the ground states is considered even more robust QGP probe as the cold 
nuclear matter effects if any cancel out and can be neglected in the ratios. The calculation 



of ratios of T states was also made in few works e.g. Il2|] in past which showed that the 
Pt dependence of such ratio would show large variations and this would be a direct probe of 
the QGP. 

In this paper, we calculate the bottomonia suppression due to color screening in an ex- 
panding QGP using the model by Chu and Matsui [isl, which takes into account the finite 
QGP lifetime and spatial extent. We start by describing the properties of quarkonia obtained 
from potential models and then give a brief description of the model which is extended to 
get the survival probabilities of T states as a function of centrality of the collisions. Finally 
we compare the model calculations with the experimental data very recently measured by 
the CMS experiment. 

II. PROPERTIES OF THE T STATES FROM POTENTIAL MODELS 

Interaction between the heavy quark and its antiquark inside the quarkonium at zero 
temperature can be described by Cornell potential jl^ which is written as 

V(r,T = 0) = ar --, (1) 
r 

where r is the separation between heavy quark pair. The parameters are string tension 



(7=0.192 GeV^, gauge coupling a=0.471, 771^=1.32 GeV/c^ and mb=4.75 GeV/c^ (l5|. The 



TABLE I. Quarkonia properties from non-relativistic potential theory [iSl . [l 



Bottonium properties 


T(16) 


, /'I D^ 
Xb(l^) 


'Y^ /o c^ 
T(26) 


1 (ob) 


Xfe(2-P) 


Mass [UeV / c J 


y.4D 


n no 

y.yy 


1 n no 
iO.Oz 


10.06 


10.26 


Radius [fm] 


0.28 


0.44 


0.56 


0.78 


0.68 


TF [fm] [15] 


0.76 


2.60 


i.y 


2.4 




Td [GeV] upper limit [18] 


2 Tc 


1.3 Tc 


1.2 Tc 


ITc 




Td [GeV] used in the present work 


1.8 Tc 


1.15 Tc 


1.1 Tc 


0.9 Tc 





solution of the Schroedinger equation for the above potential gives mass, bound state radius 
and the formation time Tp, the time needed to form a bound state after the production 
of heavy quark pairs. All parameters obtained with zero temperature potential using the 



parameter values 
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16| given above are summarized in first three rows of Table I, which 



describe well the experimentally observed quarkonium spectroscopy. 

The potential model can be extended to finite temperature with the main assumption 
that medium effects can be accounted for as a temperature-dependent potential. Instead of 
just looking at the individual bound states (at T = where quarkonium is well defined), 
one could rather obtain a unified treatment of bound states, threshold and continuum by 
determining the spectral function. Using a class of screened potentials based on lattice 
calculations of the static quark-antiquark free energy, spectral functions at finite temperature 



are calculated in a work 
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18| and it was found that all quarkonium states, except the IS 
bottomonium, dissolve in the deconfined phase at temperatures smaller than 1.5Tc. An 
upper limit on binding energy and the thermal width of different quarkonia states are then 
estimated using spectral functions in the quark-gluon plasma. Corresponding upper bounds 



on their dissociation temperatures Td 



18| are given in second last row of Table I. We used 



shghtly lower values of given in the last row to obtain a good match with measured Raa- 



III. QUARKONIA SUPPRESSION IN FINITE SIZE QGP 



The bottomonia survival probabilities due to color screening in an expanding QGP is 
estimated using a dynamical model which takes into account the finite lifetime and spatial 
extent of the system The competition between the resonance formation time tf and 

the plasma characteristics such as temperature, lifetime and spatial extent decide the pt 
dependence of the survival probabilities of T sates. We describe the essential steps used to 
develop the model which is then extended to get the survival probabilities as a function of 
centrality of the collision. 

The model assumes that quark gluon plasma is formed at some initial entropy density sq 
corresponding to initial temperature Tq at time tq which undergoes an isentropic expansion 
by Bjorken's hydrodynamics [IQi]. The plasma cools to an entropy density sd corresponding 
to the dissociation temperature To in time td which is given by 

As long as T]j/tp > 1, quarkonium formation will be suppressed. 

In finite system produced in heavy ion collision, the suppression and entropy depend on 
the size of the system. The initial entropy density is assumed to be dependent on radius R 



(decided by the centrality of the collision) of the QGP 
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as 



^ X 2\ 1/4 



so{r) = .0 ( 1 - ( ) ) , (3) 



Using Eq. (|2]) and Eq. ([3]) one can obtain the r dependence of as 



^ X 2\ 1/4 



rD(r)=rz,(0)(^l-(^-j J . (4) 

where t/5(0) is the value of To for resonances produced in the center of the system. 

Let a QQ pair is created at the position r in the transverse plane with a transverse 



momentum px and transverse energy Et = y M"^ + Pt- The upsilon formation time is Tpj 
which on equating with the screening duration To{r) given in Eq (j4]) one obtains the critical 
radius r^, which is the boundary of the suppression region as 




ro^R\l-[^^] I . (5) 

where 7 = Et/M is the Lorentz factor associated with the transverse motion of the pair. A 
bottom-quark pair can escape the screening region r^, and form T if the position at which it 
is created satisfies 

|r + ^|>rz., (6) 

where the screening region r < ro is shrinking because of the coohng of the system. Defining 
(f) to be the angle between px and r, the Eq. (Q leads to a range of (p for which the bottom- 
quark pair can escape: 

cos (p>z where z = ; — ^ , (7) 

2r{TFPT/M) ^ > 

With this we can then calculate probability for the pair created at r with transverse momen- 
tum px to survive as 



(f)(r,pT) = vr z < —1 
= cos~^z \z\ < 1 
= z>l, 

If the probability p(r) of a quark pair to be created at r which is symmetric in transverse 
plane is parameterized as 



pM=(i-(7^) 1 > (8) 



S{pr,R) = '' ^^yy;v>y,^T; _ (9) 



the survival probability of quarkonia becomes 

/(f rfr r p(r) (f){r,p T) 
IT dr r p{r) 

The survival probability as a function of centrality can be obtained by integrating over pt 
as follows 

^(iVpart) = / S{pT, R{Np.rt)) Y (pr) dpr- (10) 
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Here Y{pt) is Pt distribution (normalized to one) obtained from Pythia. The size R = 
-R(^part) as a function of centrality is obtained in terms of the radius of the Pb nucleus given 
by Ro = To A^^^ and the total number of participants A^parto = in head-on collisions as 



R{N,^rt)=RoJ^- (11) 
V jVpartO 



' part 

The initial temperature as a function of centrality is calculated by 



where Tq is the initial temperature assumed in 0-5% centrality and {dN/dt]) is the multiplicity 
as a function of number of participants measured by ALICE experiment [20]. Both ALICE 
and CMS |2l| measurements on multiplicity agree well with each other. Equation fll2p giving 
the variation of initial temperature as a function of centrality differs from the approach taken 
in the work of Ref. 22| where it is taken to vary as a third root of number of participants. 
The initial temperature Tq for 0-5% central collisions for a given initial time tq is obtained 
by 

3 3.6 (dN\ 

Here {dN/dri)o^5% = 1.5x1600 obtained from the charge particle multiplicity measured in 
Pb+Pb collisions at 2.76 TeV 2Q| and Og = 377r^/90 is the degrees of freedom we take in in 
quark gluon phase. Using Eq. f lTT]) we can obtain the transverse size of the system for 0-5% 
centrality as -Ro-5% = 0.92i?o- For = 0.1 fm/c, we obtain Tq as 0.62 GeV using Eq. (fT3|l . 
The critical temperature is taken as Tc = 0.160 GeV 

The nuclear modification factor, Raa is obtained from survival probability taking into 
account the feed-down corrections as follows, 



i?AA(3^) = ^(3^) 

i?AA(25) = /i S(25) + /2 Si3S) 

RAAilS) = gi S{IS) + S{xb{lP) + g, S{2S) + g, S{?,S)) 

(14) 
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The factors /'s and g^s are obtained from CDF measurement [23|. The values of gi, g2, gz 
and (74 are 0.509, 0.27, 0.107 and 0.113 respectively. Here it is assumed that the survival 
probabilities of T(35') and x;,(2P) are same and g^ is their combined fraction. The values of 
/i and /2 are taken as 0.50 guided by work from Ref. |24| . 
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FIG. 1. a) Measured {dN/d'r])/{Npa.j.i/2) [20'] as a function of A'^part along with the function 
{dN / dr]) / {liE?). (b) The initial temperature obtained from measured multiplicity using Eq. (|12p 



IV. RESULTS AND DISCUSSIONS 



Figured] (a) shows measured [dN / drj) / {Nps_^x_/2) [20| as a function of A^'part- The function 
{dN / dr]) / {it R"^) gives the multiplicity divided by transverse area obtained using Eq.f llip . 
Figure [T](b) gives the initial temperature obtained from measured multiplicity using Eq. (fT2|) . 
Except in peripheral collisions, the initial temperature has weak dependence on centrality of 
collisions. 

Figure [2] demonstrates the working of the model. It shows the screening radius r£, (in fm) 
in case of T(IS') and T(25') for R = 7.5 fm (corresponding to head collisions) and R = 4.08 
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(GeV/c) 



FIG. 2. The screening radius r£, (in fm) in case of T{\S) and T[2S) for R = 7.5 fm (corresponding 
to head colhsions) and R = 4.08 fm (corresponding to minimum bias colhsions) as a function of 
transverse momentum. The straight hues | r+^^-^g^ | mark the distance a bottom quark pair (created 
at r = 0) will travel before forming T(IS') or T{2S). The region below the line corresponding to 
each bound state satisfies the condition given by Eq. ([U]) and the bottom quark pair will escape the 
screening and will form the bound state. 

fm (corresponding to minimum bias collisions) as a function of transverse momentum. The 
straight lines | r+ ^-^jf^ | mark the distance a bottom quark pair (created at r = 0) will travel 
before forming T(IS') or T{2S). The region below the line corresponding to each bound state 
satisfies the condition given by Eq. ([6]) where the bottom quark pair will escape the screening 
and will form the bound state. If r is non-zero, the region where a bottomonium can escape 
screening, enlarges. 

Figure[3](a) shows the survival probability as a function of pt for T(IS'), T(2S'), T(3S') and 
Xf,(l-P) for R = 4.08 fm (corresponding to average A^part = 114 for minimum bias collisions). 
The survival probability S{pt) has a unique px dependence decided by the Tjy and Tp of each 
T state. In general, the survival probabilities of resonance states increases with increasing pt 
and become 1 at different px for different states corresponding to complete survival. Since 
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FIG. 3. (a) The survival probability as a function of pT for T(15), T(2S'), T(3S') and Xhi^P) for 
i? = 4.08 fm (corresponding to average A'part = 114 for minimum bias collisions), (b) The nuclear 
modification factor for T(IS') and T(2S') which is obtained from survival probabilities including 
feed down corrections. The solid squares are T(IS') .Raa measured in the minimum bias Pb+Pb 



collisions at ^/snn = 2.76 TeV by CMS experiment 



5i. 



T(IS') is expected to dissolve at a higher temperature it has more probabihty to survive 
the plasma region even at lower p^ as compared to the cases of other bottomonia states. 
Figure |3] (b) shows the nuclear modification factor for T(IS') and T(2S') which is obtained 
from survival probabilities including feed down corrections. The solid squares are T(IS') -Raa 
measured in the minimum bias Pb+Pb collisions at ^/snn = 2.76 TeV by CMS experiment 
The model reproduces the trend of the pt dependence of low statistics measurements of 
Raa from 2010 Pb+Pb collisions by CMS. 

Figure m^a) shows the nuclear modification factor, Raa as a function of A^part for T(IS'). 
The solid squares are measured Raa by CMS experiment in Pb+Pb collisions at ^/s^ = 2.76 
TeV {5 1 and solid circle is the minimum bias data point. Similarly, the nuclear modification 
factor, Raa for T(25') is shown in Figure |l](b). The solid line in both the figures are present 
calculations. The model correctly reproduces the measured nuclear modification factors of 
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FIG. 4. (a) The nuclear modification factor, Raa as a function of A'^part for T(IS'). The solid 
squares are measured Raa by CMS experiment in Pb+Pb collisions at -^/snn = 2.76 TeV [6] and 
solid circle is the minimum bias data point, (b) Same as a) for T(2S'). The solid line in both the 
figures represents the present model calculations. 



both T(IS') and T{2S) for all centralities using the parameters given in the Table ( HTll . 

We also calculated the ratio of Raa of T{2S) to that of T(IS') which is equivalent to the 
so called double ratio [T{2S)/T{lS)]pbPb/['^{2S)/T{lS)]pp. The double ratio has the ad- 
vantage that the effects such as initial-state nuclear effects and regeneration which we ignore 
in our calculations are supposedly canceled out. Figure |5] shows the double ratio measured 
by CMS experiment Q| along with the present calculation. The calculations reproduce the 
measured double ratio even for the most peripheral data point. 

We calculate the survival probabilities of T states and obtain the nuclear modification 
factors due to colour screening. The most important parameters in the model are formation 
time and dissociation temperatures of bottomonia states. There are reliable calculations 
of formation time obtained from zero temperature potential models which reproduce the 
bottomonia spectroscopy very well. Upper limits are available for dissociation temperatures 
which are obtained from potential models at finite temperature. We used slightly lower 
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FIG. 5. Double ratio, [T{2S)/T{lS)]phPb/[T{2S)/T{lS)]pp as a function of A^part measured by 



CMS experiment 



along with the present calculation (solid line). 



values of the dissociation temperature to get a good description of the measured nuclear 
modification factors of T(IS') and T{2S). The dynamics of the system is affected by the 
initial conditions which in present calculations are fixed using measured charged particle 
multiplicity at LHC. There can be other sources in the measured nuclear modification factor 
which we discuss here. There can be suppression due to initial nuclear suppression which 
we assume to be much smaller than that due to colour screening and hence are ignored 
in the present work. The calculations of shadowing in Pb+Pb show that it will effect the 
bottomoina yields by approximately 10 %. Thus, the dissociation temperatures obtained by 
us are still considered to be the upper limits. The bottom quark mass is 10 times higher than 
the temperatures we are considering for the system and hence the regeneration effect can be 
safely ignored in calculating nuclear modification for bottomonia. The uncertainties in the 
measurements of feed-down fractions would introduce uncertainties in the calculated nuclear 
modification factor. Finally we mention that the uncertainties arising from effects other than 
colour screening are small and supposedly will have little or no effect on the double ratio. 
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V. CONCLUSIONS 



In summary, we calculate the survival probabilities of T states and obtain the nuclear 
modification factors due to colour screening in an expanding quark gluon plasma of finite 
lifetime and size produced during Pb+Pb collisions y/s^N = 2.76 TeV. The formation time 
and dissociation temperatures of bottomonia states obtained from potential models are used 
as input parameters in the model. We used slightly lower values of the dissociation tem- 
peratures to get a good description of the measured nuclear modification factors of T(IS') 
and T(25'). The model reproduces the centrality dependence of measured nuclear modifica- 
tion factors of T(IS') and T(2S') and the double ratio very well at ^snn = 2.76 TeV. The 
trend of px dependence of low statistics measurements of nuclear modification factor from 
2010 Pb+Pb collisions of CMS are reproduced as well. The uncertainties arising from effects 
other than colour screening are assumed to be small and supposedly will have little or no 
effect on the double ratio calculations. 
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